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Abstract 



In order to explore some general features of modified theories of gravity which in- 
volve higher derivatives and spontaneous Lorentz and/or diffeomorphism symmetry 
breaking, we study the recently proposed new version of Covariant Renormalizable 
Gravity (CRG). CRG is a higher derivative theory of gravity that aims to provide a 
. generally covariant ultraviolet completion of general relativity. Its special features 

are the presence of projection operators constructed from a constrained scalar field 
■ in the action and the spontaneous breaking of Lorentz invariance, which enable the 

theory to be power-counting renormalizable. We obtain an Arnowitt-Deser-Misner 
q ' (ADM) representation of the CRG action in 4-dimensional spacetime with respect 

to a foliation of spacetime adapted to the constrained scalar field. The resulting 
action is analyzed by using Hamiltonian formalism, which was originally developed 
for constrained systems by Dirac. It is found that the theory contains two extra 
propagating degrees of freedom. One is due to the presence of second order time 
derivatives in the ADM representation of the CRG action. The other is due to 
the lack of a local Hamiltonian constraint, similarly as in the projectable version of 
Hofava-Lifshitz gravity. It is argued that CRG contains a degree of freedom that 
carries negative energy (a ghost) which will destabilize the theory. Such a pathology 
jeopardizes all higher time derivative field theories where degrees of freedom with 
positive and negative energies interact with each other, unless the given Lagrangian 
is appropriately degenerate so that there exist constraints that protect the stability. 
We conjecture that generally covariant higher derivative theories of gravity which 
involve spontaneous (constraint induced) Lorentz and/or diffeomorphism symmetry 
^ ■ breaking will in general share this problem with CRG. 



PACS: 04.50.Kd (Modified theories of gravity), 04.60.-m (Quantum gravity), 
ll.10.Ef (Lagrangian and Hamiltonian approach), 

98.80. Cq (Particle-theory and field-theory models of the early Universe) 



1 Introduction 

In the recent years modified theories of gravity have attracted a considerable amount of 
attention. These modifications of General Relativity (GR) aim to improve the behaviour 
of the theory either at high energies or at large distances. This is motivated by the fact 
that although GR is very successful in describing gravitational phenomena in intermediate 
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distances, it has become increasingly evident that it may need to be completed both at 
the high energy realm (quantum gravity and renormalizability) and at large distances 
(dark energy and dark matter) in order to achieve a more plausible theoretical framework 
and also more convincing agreement with observational data. The fact that according 
to standard cosmology only four percent of the energy content of the universe has been 
observed by means other than gravity is a huge gap in our understanding on the nature 
of gravity. 

Perturbative renormalization of GR requires us to include invariants quadratic in cur- 
vature into Lagrangian as counterterms [l|, 0|. The Riemann tensor squared term can 
be excluded in 4-dimensional spacetime due to the Gauss-Bonnet topological invariance. 
Gravitational Lagrangian consisting of the scalar curvature R, scalar curvature squared 
R 2 and Ricci tensor squared R^R^ terms is indeed renormalizable via dimensional regu- 
larization [ij . Unfortunately it includes a massive spin-2 excitation with negative energy, 
because R^ u includes second-order time derivatives of every component of the metric. 
This means the theory suffers from Ostrogradskian instability 0], an illness that hampers 
many higher derivative theories. Therefore much interest has been directed to modified 
gravity with the Lagrangian R + aR 2 p, Q , which contains only one extra scalar degree 
of freedom and it is healthy for a > 0, although the R 2 term alone is not sufficient for 
renormalizability. More generally one considers f(R) gravity whose Lagrangian is a non- 
linear function of the scalar curvature. For a recent review and references, see for example 
. These theories also have the advantage of being able to realize cosmological phases 
of accelerated expansion (inflation and present era) without additional dark components. 
Indeed f{R) gravity has an equivalent representation as a minimally coupled scalar-tensor 
theory and therefore it is practically equivalent to quintessence in this respect. Other well 
known examples of modified theories of gravity are for example Brans-Dicke theory \aL 
another scalar-tensor theory, and the relativistic tensor- vector-scalar implementation |9| 
of modified Newtonian dynamics 
Hofava-Lifshitz (HL) gravity 
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is a power-counting renormalizable field theory of 
gravity which is based on the idea that space and time scale differently at high energies, 

x-tbx, t^b z t, (1.1) 

with a dynamic critical exponent z. This enables one to modify the ultraviolet behaviour of 
the graviton propagator to |fc|~ 2z , where k is the spatial momentum. In D spatial dimen- 
sions, choosing z = D in the ultraviolet fixed point ensures the theory is power-counting 
renormalizable. Such a spacetime admits a preferred foliation into spatial hyper surf aces, 
and hence the local Lorentz invariance is broken. Essential extensions of HL gravity have 



been presented i n | 12| . |13J. Note that more general theories of the HL type have been 
proposed in Ref. [l 



Recently the so-called Covariant Renormalizable Gravity (CRG) was proposed [16 
CRG aims to provide a power-counting renormalizable field theory of gravity that is 
covariant under spacetime diffeomorphism and possesses local Lorentz invariance at the 
fundamental level. CRG aims to achieve a similar ultraviolet behaviour of the graviton 
propagator as HL gravity, but without introducing explicitly Lorentz noninvariant terms 
into the action. Lorentz invariance of the graviton propagator of CRG is, however, broken 
spontaneously at high energies. This is achieved by introducing a scalar field, which 
is coupled to spacetime in a rather complicated way, and a constraint on the scalar field 
that breaks Lorentz symmetry spontaneously. Since the CRG action contains higher-order 
derivatives it exhibits extra degrees of freedom. A new version of CRG has been proposed 
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17] . where a perturbative analysis around Minkowski spacetime showed that the extra 
degrees of freedom present in the theory do not propagate. However, we should note that 
the renormalizability of CRG, as well as of the HL theory, is assumed only based on the 
power-counting arguments. There are several potential pathologies that could ruin the 
renormalizability of this theory, such as gradient instabilities, ghosts or strong coupling. 
Since a violation of Lorentz invariance has never been observed, one could try to argue 
that CRG is a more natural modification of GR than the explicitly Lorentz noninvariant 



ones, in particular HL gravity and its generalizations (see e.g. [12|-|15]). On the other 
hand Lorentz invariance could equally well be broken explicitly at very high energies as 
long as it is somehow restored at sufficiently low energies. 

At present most modified theories of gravity should be treated as effective or phe- 
nomenological theories, since they are not derived from any compelling first principles, 
rather constructed to meet some specific purposes. Nevertheless these theories can teach 
us a great deal about the aspects of gravity, and ultimately help us in laying down the 
foundations for the next paradigm of space, time and gravity. 

Hamiltonian formalism provides a powerful tool for the analysis of constrained systems 
such as gravity. For example, it has been shown by Hamiltonian analysis that HL gravity 
can be physically consistent (N > 0) at high energies only if either the projectability 



condition is imposed on the lapse function A" [18] or the potential part of the action is 



extended with terms that contain spatial derivatives of N [19|, |20| as the vector 



a i=-Jf- ( L2 ) 

A similar result regarding the projectability of A" has also been obtained for the more 
general modified F(R) HL gravity [21]. The latter approach is the so-called consistent 



extension of HL gravity [13] , where the extra degree of freedom present in the theory [22 



has a healthy quadratic action due to extension of the action by the a, terms. In a previous 



paper [23] we studied the first version of CRG from the point of view of Hamiltonian 
analysis. It was found that there indeed exist extra degrees of freedom compared to GR, 
because there are not enough constraints to eliminate them. In this paper we analyze 
the new improved version of CRG, concentrating on the renormalizable model (z = 3) 
in 4-dimensional spacetime. In particular we are interested in whether the new model 
contains extra modes similarly as original CRG and whether the possible extra modes 
are pathological ghosts. We shall see that the results of this study will be of interest to 
other conceivable generally covariant higher derivative theories of gravity which involve 
spontaneous (constraint induced) Lorentz and/or diffeomorphism symmetry breaking. 

First the action of new CRG is introduced in Sec. [2j We obtain the Arnowitt-Deser- 
Misner (ADM) representation of the new CRG action with respect to a foliation of space- 
time adapted to the constrained scalar field in Sec. [3j The resulting action contains time 
derivative of the extrinsic curvature and powers of the extrinsic curvature up to sixth 
order. In Sec. HI introducing some additional fields enables us to obtain a first-order ac- 
tion with a kinetic part quadratic in extrinsic curvature and other first time derivatives. 
Certain solutions of the first-order action are discussed in Sec. [51 In Sec. El we analyze 
the action using Hamiltonian formalism. A few alternative sets of variables for Hamilto- 
nian formulation are considered in Sec. [71 Conclusions and some further discussion are 
presented in Sec. [HJ 
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2 Action 



We consider the new version of covariant renormalizable gravity with projectors [17]. For 
definiteness we shall consider the specific model corresponding to critical exponent z = 3 
which should be power-counting renormalizable in 4-dimensional spacetime. The action 
reads 

S 3 = J d 4 xV^g ^ - aP a »P^ fall,,,, - ^^V%V^j 



A ( -d^oy + u 



(2.1) 



where the projector is defined by 

2C/ 1 J 

Variation of the action (12. ip with respect to the Lagrange multiplier A implies the con- 
straint on the scalar field as 

^dy + Uo^O. (2.3) 
For simplicity we consider Uq to be a positive constant, as was considered in the original 



proposal [171 ]. but more generally it could be any positive function of 0. 

Choosing a solution of the equation of motion (12.31) spontaneously breaks the Lorentz 
invariance and/or the full general covariance of the action (12.11) jl71 ]. According to the 
constraint (I2.3p d^cj) is timelike and hence one can choose the direction of time to be 
parallel to <9 M at least locally. Then one finds that the modification added to Einstein- 
Hilbert action in (12. ip . i.e., the term with the coupling constant a, turns out to contain 
only spatial derivatives of the perturbation h^ v to Minkowski metric: ^'g^ = r/^ u + h^. 
In the present case (z = 3) there are six spatial derivatives of h^, while in the general case 
there are 2z spatial derivatives, yielding the desired |fc| -22 modification of the graviton 
propagator. The form of this modification to GR was specifically constructed so that time 
derivatives of h^ u cancel out in the action once the Lorentz symmetry is spontaneously 
broken and the direction of time has been fixed. This means that all the problems of higher 
time derivative theories can be avoided in this linearized formulation, quite similarly as 
in HL gravity where higher time derivatives are explicitly excluded in the definition of 
action. 

Although we concentrate on the renormalizable case z = 3 in four dimensions, actions 
corresponding to other (higher) values z could be considered in a similar fashion by using 
the results we will obtain in the following sections. 



3 Arnowitt-Deser-Misner representation 



We consider the ADM decomposition of the gravitational field [24| ; for reviews and math- 
ematical background, see 25(. Spacetime is assumed to admit a foliation into spacelike 
hypersurfaces £f of constant t. The metric of spacetime can be written as 



(4), 



(3.1) 
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where g^ v is the induced metric on the spacelike hypersurfaces H t and is the future- 
directed unit normal to S t . The ADM variables consists of the lapse function N, the shift 
vector N l and the spatial metric = 1, 2, 3). The unit normal to E t can be written 

in terms of the ADM variables as 

n, = -NV,t = (-N, 0) , n»= (n°, n l ) = -— J . (3.2) 

Thus the ADM representation for the metric of spacetime reads 

{4 W = -iV 2 + i\W% (% 0i = (% i0 = Ni , { % = g lj} (3.3) 
where iVj = gijN^ . Contravariant components of the metric of spacetime are 

(%00 L (4)„0i _ (4) iO _ ^ (4) ij _ ij _ N ' NJ ( o 4) 

y jy2 ' y y " jy2 ' y y jy2 ' ^ d I 

where g^gjk = 8 l k . The extrinsic curvature of the spatial hypersurface S 4 is defined by 

K *t = ^ - 2D « N ») ' * = 9 ijR v > ( 3 - 5 ) 

where the dot denotes the derivative with respect to time £. Quantities defined on the 
spacetime M. and associated with its metric ^g^ v are marked with the prefix ( 4 ). We 
denote the covariant derivatives on Ai and E t by V and D, respectively. For further 



details on the notation, one can also see 23 



3.1 Space-time decomposition of the action 

Instead of deriving the general Arnowitt-Deser-Misner (ADM) representation of the action 
(12. ip . as was first done in our analysis of the original formulation of CRG [23], we directly 
consider a specially chosen foliation of spacetime. The constraint (12. 3p on ensures that 
the vector d^ip is timelike everywhere, when Uq > is assumed. Therefore there exists a 
preferred foliation of spacetime into spatial hypersurfaces H t whose unit normal is given 
by 

n — / 5= — ; • (3-6) 

Then the constraint (I2.3P reduces to the condition that the normal has unit norm 

«X = - 1 > ( 3 - 7 ) 

which is presumed for the unit normal by definition. Hence we are effectively substituting 
the solution of the equation of motion ( 12. 3 p back into the action. The projector ( 12. 2 j) 
becomes the orthogonal projector onto the spatial hypersurface E t , P a fJ- = g fJ- a , which is 
defined as 

9\ = K + ^n a . (3.8) 

Now we can write 

= -^W~ n» , «9 M = -y/TXhnp . (3.9) 

From (13. 2p and (13. 9p we see that in this foliation is constant on E t , <fi = (pit). Clearly 
this construction bears some similarity to the Stiickelberg formalism used in HL gravity 



22l . l26j , where the foliation structure of spacetime is encoded into a scalar field in order to 
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achieve a manifestly covariant action and "transfer" the extra degree of freedom from the 
metric to the scalar. Here the scalar field is present from the beginning and we choose 
to work with a foliation of spacetime defined by via (13. 6p . Indeed, one could choose to 
decompose the action with respect to an arbitrary foliation of spacetime, with no relation 
to 0, but in that case the action turns out quite complicated, involving higher-order time 
derivatives up to fifth order. Note that is no longer an independent variable, but rather 
related to Uo and the lapse N. The constraint (12. 3p on reduces to 

A? 

+ U = 0. (3.10) 



2N 2 

and we can integrate it for the scalar 



</>(t) = 0(t o ) + VlU~o / dt'N{t') . (3.11) 



to 



Evidently the choice N = 1/ \/2Uq corresponds to actually choosing the scalar field as the 
time coordinate, <fi = t. The relation (I3.10P or (13. lip implies that the lapse N must be 
constant on H t too, N = N(t), because otherwise could not be constant on E t . In order 
to preserve the conditions = 0(t) and N = N(t) we restrict the symmetry under dif- 
feomorphisms of spacetime to the symmetry under foliation-preserving diffeomorphisms, 
given in the infinitesimal form as 

5t = f{t), 5x = £{t,x). (3.12) 

This is the main symmetry group of the Hofava-Lifshitz (HL) gravity. In the language of 
Hofava's theory we would say that both and N are projectable — like the lapse is in the 
projectable version of HL gravity. Here we consider (I3.12p as a partial gauge fixing of the 
diffeomorphism symmetry, which is required by the choice of working with the preferred 
foliation. 

Now the action (12 .ip can be written 

/~(4)d 
d^x^k, ^-a^(%-V n ^ + o„o„) 

x 2U (nW^V u + V^V M ) g^g"? ( { %, - V n K^ + a^a u )] , (3.13) 



where we denote V n = n^V^ and use (13.91) as V M = —y/2U n fl and apply the following 
geometrical identities 

V^n u = - , (3. 14) 

O/Lt = Vnn^ = Dp In N , (3.15) 



in order to write 



1 



2U 



<9 P 0V P V '^V u <p = - a^a v - n^V n a u . (3.16) 



The vector a M has the physical interpretation of being the acceleration of an observer 
with 4- velocity n.^. It is always orthogonal to n M , a M n M = 0, and hence tangent to S t , 
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In the action (13.131) . we recognize g^ a g v S^R^ V as the component of the Ricci tensor of 
spacetime that is tangent to E t . With the help of the Gauss relation (jA.lj) and the Ricci 
equation ( 1A.3)) it can be written 

3^7%, = Rap + KK aP - 2K all K% - -jyD a DpN + ^C Nn K aP , (3.17) 

where R a p is the Ricci tensor of the hypersurface £ t and C^n denotes the Lie derivative 
along the 4- vector Nn^ = (1, — N l ). Then we calculate the Lie derivative in (13. 17)) 

-jC Nn K a p = V n K ap + (K a » - n a a») + (k^ - n p a^ K afl . (3. 18) 

For any tensor T that is tangent to E t its Lie derivative C^ n T is also tangent to E t . This 
follows from CN n g fl a = 0. Thus we can project both sides of (13.18)) onto £ 4 and obtain 

^C Nn K aP = g^g^VnK^ + 2K m K^ . (3.19) 

Substituting ( 13 . 19)) into ( 13 . 1 7)) enables us to write 

g\g v ^% u = R aP + KK a/ 3 - ^D a D p N + g\g v ^ n K^ v . (3.20) 

Thus in the action (13. 13)) we obtain the decomposition 

g» a g% ( {A) R^ - VnK^ + a^a u ) = R a p + KK a/3 - D a a p , (3.21) 

where we can actually drop the last term in the right-hand side, —D a ap = a a a^ — 
jjD a DpN , since the condition iV = N(t) implies that a M vanishes 



a, 



N^D^N = N~ 1 g v d„N = , (3.22) 



because g° M = and d{N = 0. All the terms in the right-hand side of (13.21)) are tangent 
to E t , which is the essential geometrical implication for introducing the projectors into 
the action of new CRG. 

Using (13.21)) we can write the action (13 . 13)) as 

/r(4)p 
fayf^Vg ^ - 2aU {R af3 + KK af3 ) (nW,V u + VV„) (R^ + KK afi 

(3.23) 

What remains to be decomposed are the covariant derivatives of tensor fields that are 
tangent to the spatial hypersurfaces E t , namely 

V^V, (R^ + KK a p) . (3.24) 

In the original formulation of CRG one instead takes covariant derivatives of a scalar, 
namely the component of the Einstein tensor that is purely normal to E f , i.e. d 1 " l <fid l 'cfi^G ' ^ oc 

i^G^ = | (K 2 — + R), which are much easier to decompose. The spatial co- 

variant derivative D of a (k, /)-tensor field T on Et is given in terms of the covariant 
derivative V of spacetime as 

U V l px-pi - 9 fj,9 at 9 a k 9 pi 9 pjVa-i ft-ft ' [6. AO) 



where in the right-hand side one considers the extension of T on spacetime. This relation of 
D and V can be used in the decomposition of covariant derivatives of tensor fields tangent 
to £ t . In the Appendix IA.3| we present the calculation of space-time decomposition of 
first and second order covariant derivatives of a symmetric tensor field A a p that is tangent 
to E t . 

Here we consider the specific combination of second order covariant derivatives that 
appears in the action (13. 23ft : 

(nVV^V, + V%) A aP = g^ g^g\g\V „V p A uX - n a g^ g^n° g\V U V p A aX 

- npg u »g p ^g a a n x V v V p A aX + n a n p g u »g p ^rfn x V V V p A aX 

= D p D p A aP - K (^C Nn A a p - 2K { /A/ 3)p 

+ 2K^K p(a \A m 

+ n a (2K^D p A vP + D p K pv A% - Ka»A pP ) 
+ np (2K^D p A au + D p K py AJ - KA ap a») 
+ n a n^2K^K/A up . 

(3.26) 

When we multiply (I3.26P with another symmetric tensor field B a/3 that is tangent to £ t , 
we obtain 

B a ? (nW p V u + V P V P ) A a p = B ij [D k D k A ii - K^-C Nn Aij 

V N (3.27) 

+ 2KK i k A kj + 2K ik K kl A lj ) . 



In the right-hand side of f!3.27j) all the tensor fields are tangent to S t , which enabled us 
to write the contractions over spatial components. Substituting 

Aij = Rij + KKij , B ij = R ij + KK ij (3.28) 

into (I3.27P gives the ADM representation of the action (13.131) 



D k D k (R^ + KKij) 

. (3.29) 



/( K — K 2 -\- R 
d'x^gN 2 f +K - 2aU + KK ij ) 

- K^C Nn ^ + KKij) + 2KK k (R kj + KK kj ) + 2K ik K kl (R l3 + KK b 

Here we have also included the usual decomposition of the scalar curvature of spacetime 
(1A.4|) and dropped the divergence terms; we assume the spatial manifold is compact and 
that it has no boundary. The first thing to note is that the Lagrangian of the new CRG 
given in (13.291) contains kinetic terms involving the extrinsic curvature to the sixth power, 
while conventional gravitational Lagrangians only include kinetic terms quadratic in the 
extrinsic curvature. The action also contains second order time derivatives of the metric in 
the term involving the Lie derivative of the extrinsic curvature. Namely in the Lagrangian 
we have a term that contains 

K (R ij + KK ij ) C Nn (Rij + KK^) = K (R ij + KK lJ ) [d t (R t] + KK t] ) 

- C N (Rij + KKij)} , { ' ' 
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where denotes the Lie derivative along the shift vector N l in £ t . Space-time decom- 
position of the Lie derivative C^n of a covariant tensor tangent to £ t can be found in 
Appendix IA.3I 



Recall that in the linearized treatment [17| the part of the action with the coupling 



constant a contained only spatial derivatives of the perturbation to Minkowski metric. 
Here we obtain that our ADM representation of the same part of the action contains first 
and second order time derivatives. Therefore already at this point we can expect to obtain 



some results that will differ from [17] , in particular due to the presence of the higher time 
derivative term (I3.30p . 

ADM representations of CRG actions corresponding to values of the critical exponent 
z other than z = 3 can be written by using the results (13.211) and (I3.26p . It is found 
that in general such actions contain kinetic terms involving the extrinsic curvature and 
possibly its time derivatives to the 2z-th power. Higher-order time derivatives are present 
because every instance of the derivative operator decomposed in (I3.26P adds one more time 
derivative. The only exception is the nonrenormalizable case z = 2, where the action is 
like (13.231) but without the derivative operator decomposed in (I3.26p . which involves only 
first order time derivatives albeit its kinetic part is still modified substantially. We will 
not present the ADM representations of the actions for higher z here since generalization 
of the present results is straightforward. 

ADM representations of other conceivable generally covariant higher derivative theo- 
ries of gravity which aim to be power-counting renormalizable by involving spontaneous 
(constraint induced) Lorentz and/or diffeomorphism symmetry breaking will likely share 
some characteristics with CRG. In particular, it is unlikely that all higher time derivatives 
could be cancelled by such construction. However, constraints that also contain higher 
time derivatives might just be able to serve such purpose. 



4 First-order Lagrangian for Hamiltonian formalism 



A Hamiltonian formulation for higher derivative theories with regular Lagrangians was 
first developed by Ostrogradski [271. Some decades after Dirac developed Hamiltonian 
formalism for constrained systems |28|, |29| it was generalized to higher derivative theo- 



ries [3CH32) . Soon Hamiltonian formulations of higher derivative theories of gravity were 
constructed for the first times 33M35I]. Hamiltonian formulation of actions that involve 
higher-order time derivatives requires one to introduce a pair of new independent variables 
for each higher-order time derivative of a variable. There generally exist several different 
choices of such additional variables, which each yield a different Hamiltonian formulation 
of a given higher-order action. Such Hamiltonian formulations are connected by canon- 
ical transformations 35) and hence they are classically equivalent. But those canonical 
transformation can be highly nonlinear. Thus there is no guarantee that the Hamiltonian 
formulations remain equivalent after quantization. 

Regular Lagrangians that depend on higher-order time derivatives of dynamical vari- 
ables are known to possess degrees of freedom which carry both negative and positive 
energies 



27 



In an interacting higher derivative field theory such ghosts necessarily 
destabilize the theory when the Lagrangian is regular (nondegenerate), since any state of 
the system can and will further decay into excitations with compensating negative and 
positive energies; for an example, see 36]. Higher derivative theories whose Lagrangian 



are singular (degenerate) can sometimes avoid the Ostrogradskian instability. Theories 
which possesses continuous symmetries are always degenerate, in particular gauge theories 



9 



such as gravity, and hence they have a chance to avoid the instability. Thus for higher 
derivative theories of gravity, the existence and behaviour of ghosts has to be checked 
in each theory. The same applies to other pathologies such as strong coupling of extra 
degrees of freedom. 

A relatively simple example of a higher derivative theory of gravity is provided by 
f(^R) gravity, whose Lagrangian is an arbitrary function of the scalar curvature of 
spacetime. The scalar curvature is second order in time derivatives (IA.4j) . As a result 



the theory contains an extra degree of freedom compared to GR, whose Einstein-Hilbert 
Lagrangian is just ^R. For a recent review of f(^R) gravity, see for example [7]. f(^R) 
gravity has several known Hamiltonian formulations based on different choices for the 
higher-order variable associated with the extra degree of freedom. One can regard the 



scalar curvature of spacetime as the additional dynamical variable [341 ] . Since the second 



order time derivative in the scalar curvature (1A.4|) appears as K, we can follow the ap- 



proach of |35| by regarding the trace of the extrinsic curvature as an alternative scalar 
variable. A third alternative is provided by the fact that f(^R) gravity is equivalent 
to a scalar-tensor theory where the scalar is minimally coupled to GR. The extra scalar 
degree of freedom in f(^R) gravity is indeed stable if the potential is well behaved. For a 
discussion of the Ostrogradskian instability and the reason why f(^R) gravity is able to 



avoid it, see I4fl. For a recent review of the aforementioned three Hamiltonian formulations 



of f(^R) gravity, see [37 



In the action ( I3.29p . there is a term (I3.30p that contains a time derivative of each 
component of the extrinsic curvature K%j. That is a second order time derivative of each 
component of the metric. Therefore we shall introduce two new symmetric rank 2 tensor 
variables. The way in which these new variables are defined is a matter choice. Since 
our action contains the second order time derivatives in K and it would be natural 



to follow the approach originated in [35| where are taken as additional independent 
variables. We, however, choose alternative variables that simplify the action further. 
In order to obtain an action whose Lagrangian is quadratic in the extrinsic curvature 
and the first time derivatives of the additional higher-order variables we shall introduce 
an additional independent symmetric tensor variable Qj, which is related to the metric 
variables by 

dj = Rij + KKij (4.1) 

This choice also enables us to get rid of terms that involve the time derivative Rij of the 
spatial Ricci tensor. In order to enforce the relation (14. ip another symmetric tensor field 
will be introduced as a Lagrange multiplier. We replace the action (I3.29p by 

S 3 = J d A x^fgN KtjRt3 ~f 2 + R - 2aU Q g ik gi l ( kl (o k D k ( tJ - K^C Nn Q 3 

+ 2KK k Q kj + 2K lk K kl C h ) - o£* - R tJ - KK tJ )j . (4.2) 

Variation of the action with respect to £* J gives the equation (14. ip . Substituting it back 
into the action would give the original action (I3.29p . Since the action does not contain any 
time derivatives of £ lJ these variables have an auxiliary character; they are not expected 
to carry any physical degrees of freedom. In the action (14.21) . the term that contains a 
time derivative of Qj can be written as 

2aU Kg lk g> l ( kl £ Nn Q J = 2aU Kg* k gJ l ( kl - C N Cii) , (4-3) 
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where Cn denotes the Lie derivative along the shift vector N l in £ t . Thus at least some 
of the variables Qj are expected to carry propagating degrees of freedom. 

In the actions f !3.29|) and (14. 2 j) A" had to be restricted to a function of time only as was 
discussed in Sec. [3j In these actions, it would be possible to promote the lapse function 
to possess dependence on both space and time, since the symmetry under A3. 1 2|> would 
be retained. But the resulting action might possess very different properties compared to 
the ADM representation we have here derived from the CRG action (12. ip . Hamiltonian 
structure of such a generalized action would certainly bare some differences, in particular 
the existence of a local Hamiltonian constraint, and its study would be interesting. Here we 
shall, however, constrain ourselves to the study of the action ( 14. 2 p with the projectability 
of A" implied by the constraint on the scalar field ( I3.10p . 

Two alternative sets of variables for Hamiltonian formulation of new CRG are discussed 
in Sec. [7J 



5 On solutions for the first-order Lagrangian 

Let us consider solutions for the first-order Lagrangian, i.e., configurations for which the 
action (14. 2 p is extremal. When the Lagrange multiplier fields £ y ' are set to zero, the 
equations of motion obtained by varying the higher-order variables Qj are linear and 
homogeneous in Qj. They have a static solution Qj = provided that the initial value for 
dj is zero everywhere on the t — hypersurface So- Then the equations of motion ( 14. ip 
obtained by varying £ y give 

Rij + KKij^O, (5.1) 

which imply that the scalar intrinsic and extrinsic curvature are related by R = —K 2 . 
Substituting these solutions back into the action gives 

/K-KV —?K 2 
<*V^ 2/g2 2K (5.2) 

with the intrinsic and extrinsic curvature related by (15. ip . This kind of action with no 
curvature dependent terms in the potential is called ultralocal gravity. Such an ultralocal 
theory of gravity where the kinetic part of the action is identical to that of GR and the 



potential consists only of the cosmological constant was originally proposed in [38] . The 
Hamiltonian of ultralocal gravity has the same form as in GR, written in terms of the 
ADM variables as 

H = J d 3 x (NU + WUi) , (5.3) 

where Hq is the Hamiltonian constraint and Hi is the momentum constraint. The ap- 
pealing characteristic of ultralocal gravity is that the algebra of constraints is a true Lie 
algebra, because the Poisson bracket of Hamiltonian constraints 7io(x) and l-Lo(y) is zero. 
In GR the Poisson bracket of Hamiltonian constraints is equal to a linear combination 
of momentum constraints Hi(x) with a field-dependent "structure constraint". Thus the 
constraints of GR lack Lie algebra structure. The momentum constraint of ultralocal is 
identical to that of GR and hence the rest of the constraint algebra of ultralocal gravity 
coincide with GR. 

Note that the kinetic part of the action (15.21) is not identical to GR since its second 
term is modified by the factor 2. Rather this action formally corresponds to the ultralocal 



case of HL gravity [ll[ in the limit z — > 0, A — » 2 and with vanishing cosmological constant. 
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Here A is the coupling modifying the kinetic part of the Lagrangian: K^K 1 ^ — XK 2 . The 
only difference is that here the spatial Ricci tensor is given in terms of extrinsic curvature 
by (15. ip . A special case of this kind of solution is static flat spacetime, so that Ry = 
and Kij = satisfy (15. ip . Thus we confirm that CRG action has a static flat solution. 
However, this does not necessarily mean that it is a stable vacuum state. 

When the initial geometry of spacetime deviates from (15 .ip even slightly, the dynamics 
of spacetime produced by the action (14. 2 p is very different compared to the case discussed 
above. The equations of motion obtained by varying and Qj are highly nontrivial. 
Indeed the dynamics of the metric and the higher-order variables Qj is very complicated. 
Thus we do not attempt to solve the equations of motion in the general case. Instead we 
seek to understand the system through Hamiltonian analysis. 



6 Hamiltonian analysis 

Let us analyze the new CRG action (14. 2p by using the Hamiltonian formalism generalized 



for constrained systems |28l . l29f . For reviews, see [39] . 

In Sec. El we saw that when Qj = 0, i.e., when ( 15. ip holds, the action reduces to an 
ultralocal special case of HL gravity ( 15. 2p . Now we aim to understand the structure and 
dynamics of the theory when Qj is nonzero. We assume that at least two components of 
Qj are nonzero, while the rest of the components can attain any values. In principle it 
makes no difference which of the components are chosen to be nonzero. We choose one of 
the nonvanishing components to be the trace of Qj due to notational elegance. 

First we define the canonical momenta. Since the action (14 .2p is independent of the 
time derivatives of N, N l and their canonically conjugated momenta, p^, Pi and p^, 
respectively, are primary constraints: 

pi(x)&0, p^(x)^0. (6.1) 
The momenta canonically conjugate to gij and Qj are defined by 

'K'J - g ij K aU 



ij ^3 i- 

P = \/g 

09ij 



+ -jj-9 [Cm - L> N Q k i )g g Cn 



2k 2 

- 2aU (g ll K kl + g ik K» + g^K* + g lk tf [ K) C, km Q n g r 



+ ^(g l >e i K kl + e 3 K)\, (6.2) 
P? = = V9^U Kg* k gi l ( kl . (6 . 3) 

We shall adopt a convention where the trace component of a tensor or a tensor density is 
denoted without indices and the traceless component is denoted with the bar accent. For 
instance according to this convention we decompose the variables Qj and p 1 ^ as 

Qj = Qj + ^gijC , Vl = Vl + ^9ijPc , (6.4) 

where Qj and p 1 ^ are the traceless components and ( = g^Qj and p^ = g%jP % ^ are the 
trace components. The same notation can be used for the other variables and tensors in 
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general. Taking the traces of the momenta (16.21) and (16.31) gives 
K 3aU n 



P = V9 



- 2 + - CnSv) 9 tk 9%i ~ aU (lOK* + 2Kg*) ( lk Q l9 kl 



+ ?L(3e j K ij +£K)\, (6.5) 
p c = ^2uU K(. (6.6) 

where we use the aforementioned convention to denote p = gijP^ , £ = etc. Assuming 
has a nonvanishing trace, £ 7^ 0, and the coupling constants a ^ and Uq > 0, we 
solve (16. 6p for the trace of the extrinsic curvature 

K=±^ (6.7) 

Additional primary constraints can be obtained by substituting (16. 7p back into (I6.3p . The 
trace of the resulting equation is a trivial identity, but its traceless part expresses p]? in 
terms of other variables. Thus we obtain five new primary constraints 

ik jlJ P£ 



n v =p?-sV0BT«o I (6.8) 



which form a symmetric traceless tensor density. This primary constraint tells us that 
only the trace component of the momenta p 1 ^ is an independent variable. In order to solve 
the traceless component Kij of the extrinsic curvature we separate it from the already 
solved trace component (I6.7P as 

A-,=A-, + i 9 ,A' = A-, + -L^|. (6.9) 

Then we substitute (16 .7p and (I6.9P into (16. 2p and (I6.5p . The resulting equations can be 
solved for K^. We obtain the equation 

1 5 ( ik jl _ 1 ij kl\ /- A „mnPC ^ ( cij „ij £\ P< 



\^9 r ,i / /,■/ 
2k 2 ' 



;F^K kl , (6.10) 

where we have defined 



F ijkl = 1 (gikgjl + g il gj kj _ gK 2 aUo ^gm(i g m g l)n _ ^g^^ . (6.1l) 

pijki - g symmetric j n it s f| rs t two indices and also in its last two indices: 

pijkl _ pjikl pijkl _ pijlk 

F l i also has the properties 

gijF ijkl — g kl , (6.13) 

g kl F^ kl = g* - 8K 2 aU U k g^ 1 - \<P<1*\ CUmg mn , (6. 14) 
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where the first result is directly related to the fact that both sides of (I6.10p are traceless. 
Finding the inverse F^ kl of F tjkl , such that F^ kl F Umn = 5^8$ = 5^5^ , enables us to 
solve Kij from (I6.10p . 

In order to prove that the inverse F^ kl of ( 16. lip actually exists, and consequently that 
no more primary constraints are required, we shall construct F^ kl . First note that F~- kl 
must have similar symmetries ( I6.12p as F^ kl . We can define F^ kl as a power series in the 
coupling 2K, 2 aU with the terms depending on the spatial metric and the field Qf 

1 oo 

FA = ~ (3*9* + gugjk) + £ (2« 2 «^o) 1 c$ mini - mjBJ 

1 j=i (6-15) 

v t t n° lPl ■ • ■ C C n° jPj 

where the coefficients (j^^ linv " rnjnj depend only on the spatial metric and possess the 
symmetries of F^ kl . These coefficients are defined as follows. The first coefficient is 

eg™ = 4 (s$g j){k S% - \ gij 5p^j (6.16) 
and the rest of the coefficients are defined recursively by the formula 

rtWrnmi-mjn; _ n(J-l)mim-mj-inj_i / ( „mj(oep) rnj „°PxmJX n J] (P. 

^ijki —^ijop 4 iy °(k°i) ~ g~y % °i) ) k®- 1 - 1 ) 

for every order J > 1. The coefficient Q < ^^ aiqi '" mjqj Q f order J in (I6.17P vanishes if the 
coefficient Cf- kl 1 ^ mi?1 m J-i<?./-i Q f orc i er J — 1 i s proportional to gu- In that case the series 
would be truncated. However, it seems that such a coefficient does not appear and hence 
the recursion (16.171) will proceed ad infinitum for coefficients with increasingly complicated 
index configurations. For example the second order coefficient is 



^,(2)m 1 q 1 m 2 g2 
^ijkl 



n I n m i m 2 Xli r<?2 _i_ _gim2 rmi r92 _ " Jmxqx jrm 2 r<?2 

gi/y 1 y °(fe°z) "t- y °(fe °«) gy °(fc °j) 



(6.18) 



g vC$? ini '" mjnj =0> ( 6 - 19 ) 



Note that for every order J 
which implies 

^>«h = 0w • (6-20) 
Thus the traceless component of can be written in terms of the canonical variables 

as 

^ = ^iP«, (6.21) 

where we introduce 

^ H = P ij - \g lJ P + 1 (y V - \sf*s/*) CUmg mn ^ - ^- (V - \g ij i) ^ 

I.I u / /-k f\m 1 JW/- Amn\ ^ tW^C 



— «i I _ I /-re Aim _n i V A mrt I £_i. __£'' c '£_2. 
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for the purpose of shortening the following expressions for the Hamiltonian and its con- 
straints, which are quite complicated. From now on we may raise and lower spatial indices 
with the spatial metric, e.g., (\ = 9^(jk an d J — 9 %k 9^Ckh but always keeping in mind 
that the metric has to be written explicitly when Poisson bracket is evaluated, similarly 
as it has to be done with the trace components, e.g., ( = g^Qj- Together with (I6.9p . 
f)6.2ip and (16. 22p express the extrinsic curvature (13 .5p in terms of the canonical variables. 

Then we shall introduce the Hamiltonian. In the Lagrangian density C of the action 
(14. 2 p the term that contains the time derivative of Qj can be written 



2aU 0y /gg ik g jl ( kl K£ Nn ( ij = p^ (4/ ~ £*r&j) > 



(6.23) 



where the first term simply cancels p^Qj m the Legendre transformation from Lagrangian 
formalism to Hamiltonian formalism. For the metric we obtain 



N 



N 



P P( 



y/g3aU ( 



+ 2p i W {i N j) 



We define the total Hamiltonian as 



(6.24) 



H= I d 6 x [ P v 9ij + /A',. + X NPN + + \V p t. + v l3 W - C 



d 6 x [NU + N l Ui + X nPn + X'pi + Af4- + VijU 



V-'i 



(6.25) 



where A at, A j , A^ and Vij are Lagrange multipliers and we have defined 



Ho — —p 
V9 



+ 
+ 



P Pc 
3aU ( 



l2n 2 a 2 U$ V C 
20k 2 



v 4\ + 



9^.2 rp-1 pkl im jnrp-1 pop 
^ r ijkl r tJ tJ r mnop r 

2 

PC 

C 



j iklm 



c 



nopq 



PC 



(6.26) 



and 



Hi = -2 9lJ D k p> k + diQkpf - 2dj 



^9ijd k p jk - {2djg ik - <% jfc ) p> k - 2(^8^ - {2d£ ik - di( jk ) p J ( 



(6.27) 



Observe that the Hamiltonian only depends on the traceless components p l j! of the mo- 
menta canonically conjugate to Qj through the primary constraint (16. 8p . which also means 
that p 1 ^ does not appear in any other constraint. 

The primary constraints (16. ip and (16. 8 j) must be preserved under time evolution gen- 
erated by the total Hamiltonian (I6.25p . For this purpose we introduce the secondary 
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constraints 

J3 



$o = J d A xH « 0, (6.28) 

n t (x)^0, (6.29) 

*y(x) = {pf i (x),$ o }«0,, (6.30) 

5&(x) = {fi^aO,$ }«0, (6.31) 

which ensure that the primary constraints pn, Pi, p\j and IP- 7 are preserved in time, 
respectively. The constraint ( 16 .30 j) can be explicitly defined as the following symmetric 
tensor density 



(6.32) 



where we denote 



J.. _ JIp-l _ p-l km lnjp-1 pop , 5 ^fc^ m p-1 PC 
Yi] — P r klij 1 klijif y r mnop r ^ l r kmij £ 

The trace of the constraint \Py is independent of £ tJ '. Indeed it can be written as 



(6.33) 



This implies that 



* = ^(t) (6-34) 
{n*)Mv)} = {*«(»),PL(V)} =0. (6.35) 



The traceless component of the constraint depends on the traceless component of 
but it is independent of the trace component £. Therefore we shall also decompose each 
of the constraints p\ • and into a traceless component and a trace component: 



/';, = Pij + o9ijP* 



3^ 

^ = ^ + . (6.36) 

The explicit form of the constraint rP/ 7 defined in (16.3ip is a very complicated symmetric 
traceless tensor density. We shall not write it down here since its explicit form will not 
be used in our analysis. Instead when we evaluate Poisson brackets between IVjj and the 
other constraints we shall take advantage of its definition as the Poisson bracket of PP- 7 
and $o- 

Then we must ensure that every secondary constraint is preserved in time. First we 
shall consider the constraints "Hj. We introduce a global smeared version of the constraint 
"Hi as 

*s(x*) = J d^Ui, (6.37) 
16 



where X % {i = 1> 2, 3) are arbitrary functions on £ t which vanish rapidly enough at infinity. 
The Poisson brackets of the constraint (I6.37P with the canonical variables are 

{$s(x k ),9ij} = -X k d k 9ij ~ d l x k gkj ~ djX k gki = C x 9ij , 
{®s(x k ),P i3 } = -d kX k p ij - X%p ij + d kX l p kj + d kX j p ki = C x p ij , 
{Mx% &j} = -X k d k (ij - d iX k Ckj ~ d jX k (ki = C x (ij , 

s{x k )Jl) = ~dkX k p l l ~ X k d k pl + d kX *pf + d kX ^pf = C xP l l . 

The constraints (16.371) also satisfy the diffeomorphism algebra 

{<S>s(x l ),$sW)} =^>s(x J d jV t -V 1 d,x 1 )- (6.38) 

Thus we identify (16.371) as the momentum constraint that generates diffeomorphisms in 
the spatial hypersurface H t for the dynamical variables g^-, p % \ Qj and p^ . In fact we 
can extend the momentum constraint (I6.37P to a full generator of spatial diffeomorphisms 
with the help of the primary constraints (I6.ip . We redefine 



d 3 x X %i , 



(6.39) 



where 

$ l ='H l + C NPl + 2d jP ie k + P%d£ k + 2p\ k d ] e . (6.40) 
Note that the generator for N and pn vanishes, 

= J = 

since these variables are spatial constants. The extension (16.391) of the momentum con- 
straint simply amounts to rewriting the Lagrange multipliers of the primary constraints 
flSTTj) as 

A i = v* + C N N { , A| j = vf + C N C j , J d 3 x\ N = v N , (6.41) 

where v\ Vg and are arbitrary. The constraint (I6.39P evidently generates time- 
dependent spatial diffeomorphisms for all the variables. Therefore we obtain 

{Mx%n } = -x'diHo - d iX l H Q . (6.42) 

In other words, "Ho is a scalar density on the spatial hypersurface £j. We again emphasize 
that the lapse N depends only on time. As a result there is no local Hamiltonian constraint 
Ho but only the global one $o defined in (16.281) . As a result we clearly have 

{$ ,$ } = 0, (6.43) 
{$5(^)^0} =0. (6.44) 

It is now evident that the momentum constraint (I6.40p is preserved in time, since (16.431) 
and (16.441) and all the rest of the constraints in the Hamiltonian transform as scalar or 
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tensor densities under spatial diffeomorphism. Indeed we can write the Hamiltonian as a 
sum of the constraints 



H = N$ + $ S (JV*) + v nPn + J d 3 x (v% + jip% + Vijtt*} . (6.45) 

We could also include the 11 secondary constraints and ITjj into the Hamiltonian with 
arbitrary Lagrange multipliers. However, those 11 Lagrange multipliers would already be 
fixed by the 11 consistency conditions for the primary constraints and ft 1 - 7 . Thus the 
inclusion of those constraints is of no benefit to us. Hence we leave them out from the 
Hamiltonian. 

The global Hamiltonian constraint $o is preserved in time due to the constraints ■, 



{$o, H} « - y d 3 * (r-M', ; + ^fl^) « . (6.46) 

Next we have to ensure that the secondary constraint is preserved in time. The 
Poisson brackets j 1 !^.,, if} and must both be either zero or a constraint. The 

consistency conditions can be written as 

N{%(x),<S> }+ J d 3 y (yl\y) [^(x)^^)} + v kl (y) {%(x), Ii kl (y)}) = (6.47) 
and 

N $ } + 3a^g^^(x) = (6.48) 

for the traceless component ^ and the trace component respectively. In (16.481) . we 
used (16.351) and 

fl^(y)} = 3a^^-(x)5(x - y) . (6.49) 

The consistency condition (16.481) can be solved for one of the 5 independent components in 
the Lagrange multiplier v^. Let us denote the specific solution to (16.481) by v'^, which still 
contains 4 independent arbitrary components. The solved component of v'^ is proportional 
to N. The homogeneous part of (16.481) has only the nontrivial solution. Consistency of ^ is 
assured by substituting the solution = v'^ into the Hamiltonian. Then the consistency 
of the constraint ^ is assured by solving the inhomogeneous linear equation (I6.47P for 
the Lagrange multiplier . None of the three Poisson brackets in (I6.47P vanishes, not 
even weakly. Explicitly the second Poisson bracket can be written as 



+ F£ - lg kl F^ mn 9 mn ) (^) (x)S(x - y) 



(6.50) 



where we denote 



A — _ n mo n n P ( — —n P" 1 rfl T I 

s~*-ijkl klij r mnijy i/ I r opkl ^kl r opqriJ I 

+ %K 2 aU Q C mn F-l ij e n g oq (f^ m - \gkiF-U 



(6.51) 
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Thus we obtain 

g kl (y){%(x),pl l (y)} = 0, (6.52) 
where we also used the property (I6.20p . Clearly the homogeneous part of the equation 

dsazD 

J d 3 yv kl (y) {*v(*).Pk(»)} = ( 6 - 53 ) 



vf = \v^g kl , (6.54) 



has the solution 

where the trace component v% of the Lagrange multiplier vf is still an arbitrary field. This 
is the only nontrivial solution to the homogeneous equation (16.531) . The specific solution 
to the inhomogeneous equation (I6.47P has the form 

vf = Nwf[ gij J j Xii,Vt,C j AjlN] , (6.55) 

where wf is a traceless functional of the listed variables and of the Lagrange multiplier 
v[j divided by N. Formally it can be written as 

w*(x) = - [ d 3 yB^ kl (x,y) {*«(j/),$ } 

, (6.56) 
- J J d 3 yd 3 z^^-B i ^ n (x,y){^ mn (y),tl kl (z)} , 

where B^ kl (x,y) is the inverse to (I6.50p . i.e., it satisfies 

J d 3 yB^ l (x,y){^ kl (y),pi n (z)}=5^5(x-z). (6.57) 

The general solution to (I6.47P is the sum of the specific solution (16.551) and the solution 
(I6.54p to the homogeneous equation (I6.53[) : 

vf = Nwf [ 9ii J\ C«,PC v[ 3 /N] + l -v i9 kl . (6.58) 

Inserting this solution into the Hamiltonian (I6.45P ensures the consistency of ^>ij. The 
Hamiltonian is now written as 

H = N$ + N J (Pxw^ + <5> S (N' 1 ) + v nPn + J d 3 x (v^ + ^v^ + v'^U^ . (6.59) 

Lagrange multipliers that are arbitrary or at least contain some arbitrary components 
are denoted by v and the multipliers denoted by w have been solved entirely in order to 
ensure the consistency of constraints, e.g., in (I6.59P w 1 ^ is given by the specific solution 
(I6.56P to (I6.47P and v'^ is the solution to (I6.48j) with four arbitrary components left. 

As a last step in Dirac's algorithm we have to ensure the consistency of the secondary 
constraints (I6.3ip in time. We require that 



{^(x^HjnNl&^x^oj + N J d 3 ywf{y) {tV^x),^)} 

+ / d 3 y l -v^y){^ I (x)^(y)}+ [ d 3 yv' kl {y) {^(x), U kl (y)} (6.60) 



19 



must be either zero or a constraint. In this equation is given by (I6.56p . For the calcu- 
lation of the Poisson brackets of tljj with and p^ kl we can use the definition of fP^cc) as 
the Poisson bracket {n iJ '(aj), $o}) an d the Jacobi identity for the Poisson bracket. First 
we obtain 

{fiftC^.piCy)} = -{{$o, P l l (y)}^(x)} - *o} 

= {* k i(y),fP j (x)} (6.61) 
= {* fc Ky),n < ^(aj)} + i(/ w (y){*(y),n«(aj)} , 

where we also used the independence of IP 3 ' on For the trace component p^ we obtain 

{n? J (*) ) ^(»)} = -{{* y(w)},fl°>)}-{{p f (w),^(*)},*o} 

= - {{$ ,<7 fc, (v)} ,n«( a j)}pl I (y) + {*(i/),n«(aO} (6.62) 

w {*(!/), n«(x)} , 

where the last expression is given in (16.491) . For the traceless component pf. we obtain 

{fi&aO.j&G/)} = {n£0«0.J&(v)} - ^ {n?i(*),«w(y)}^(y) 

-^(i/){n&(*),^(y)} (6 ' 63) 

which is a complicated nonvanishing expression. The first and the last Poisson brackets 
in (I6.60p are neither zero nor a combination of the constraints. Indeed they are very 
complicated expressions. Thus we see that the 5 consistency conditions for Wjj should be 
solved for the Lagrange multiplier v% and the 4 arbitrary Lagrange multipliers left in v' kl . 
We obtain the following equations for them 

N{^j(x),^ } + N J J d 3 yd 3 z{n ij (x),^ kl (y)}B klmn (y,z){^ mn (z),^ } 

+ ay /g?L v £ x ) + J d 3 yv' kl (y) {fl«(a!),n"(y)} = 0. (6.64) 

The general solution to (I6.64p has the form 

v^ = Nw[g ij}P ij } Cij,PoC j ] , (6g5) 

v'h = N ™ki [ga , v 13 , dj ,Pc,C J ] + fhi [gn , v %3 , Cy , Pc > C J ] , 

where w and w k i comprise the specific solution to (I6.64p (divided by N), h k i is a possible 
nontrivial solutions to the part of (I6.64p that is homogeneous in v' kl , and / is an arbi- 
trary functions of time. Because of the very complicated form of the equations (I6.64p . 
we have not been able to rigorously establish the existence of neither the specific solu- 
tion nor the homogeneous solutions. Still we assume that a specific solution exists. It is 
unclear whether nontrivial solution h k i to the homogeneous part of (16.641) exists. Con- 
ceivably there could even be several solutions to the homogeneous part. The existence 
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of a solution to the homogeneous equation, would imply the existence of an extra gauge 
symmetry associated with a global (integrated) first-class constraint, a generator of the 
gauge transformation. We, however, suspect that such a gauge symmetry does not exist 
in the theory. Nevertheless if homogeneous solutions do exist they can be dealt with by 
introducing some global gauge fixing conditions. Hence we shall continue our analysis as 
if no nontrivial solution to the part of (16. 64p that is homogeneous in v' kl exists, i.e., we 
assume hu — 0. 

The Hamiltonian that ensures the consistency of all the constraints in time is 
H = N% + N j d?x (w'lplj + ^ V + w^j + $ 5 (AH + v nPn + J d 3 xv l Pl , (6.66) 

where Nw^ is the (traceless) specific solution f!6.56j) to the inhomogeneous equation (I6.47P 
and Nw and NiBij are the specific solution to the inhomogeneous equations (I6.48P and 
(I6.64p . Since Hamiltonian is always a first-class quantity and (16.66P is a sum of constraints, 
we can conclude that the Hamiltonian is a sum of first-class constraints. That is exactly 
what one expects from a system with time reparameterization invariance. We already 
knew that the momentum constraint Hi, or rather its generalization $j in (I6.40p . and the 
primary constraints p^ and pi are first-class constraints. Now we can see that the linear 
combination 

$ + J d 3 x + ^wpt + Wij&n = J d 3 x (n + //f/c ; + ^w P t + wyirn (6.67) 

is a first-class constraint. All these first-class constraints are associated with the symme- 
try under foliation preserving diffeomorphism: invariance under time-dependent spatial 
diffeomorphism and time reparameterization. In addition a global first-class constraint 
associated with each solution to the homogeneous part of the equation (I6.64p could con- 
ceivably exist, although we suspected that such solutions do not exist. 

The nature of the rest of the secondary constraints remains to be elaborated. Namely 
we are interested in whether some linear combination of the secondary constraints is a first- 
class constraint. That is in addition to the known first-class combination (I6.67p . which 
is indeed the only one that can include $o- Under the Hamiltonian ( 16.45)) we remarked 
that including ^ and H^ in the Hamiltonian with Lagrange multipliers would result to 
the fixation of the Lagrange multipliers due to the consistency conditions of the primary 
constraint pf ■ and This shows that every linear combination of these constraints is a 
second-class constraint. Consequently there is no need to extend the Hamiltonian (I6.66P 
with further secondary constraints. 

Finally we can seek to identify the physical degrees of freedom. First let us count 
the number of propagating physical degrees of freedom (physical d.o.f.) by using Dirac's 
formula: 

^(physical d.o.f.) = - [^(canonical variables) — 2 x ^(first-class constraints) 

— # (second-class constraints)] . (6.68) 
We have 42 as-dependent canonical variables (N\ gij, Qj, £ y and their conjugated mo- 

nf ] I ' ' \\l 



menta), 6 first-class constraints (pi, $«) and 22 second-class constraints (pf-, n jJ , \I>i 



42 — 12 — 22 

#(physical d.o.f.) = = 4 . (6.69) 
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That is 2 more physical degrees of freedom than in GR. For ^-independent zero modes 
there exist two more variables N,p N and two more first-class constraints pjq and (I6.67p . 
thus yielding 3 physical degrees of freedom, one more than in GR. Let us further seek to 
understand the nature of the extra modes in the theory. 

The constraints of the theory enable us to regard some of the canonical variables 
as being dependent on other variables. By analysing the dependencies of the variables 
we seek to identify the physical degrees of freedom. First let us consider the second- 
class constraints. We assume that any conceivable first-class constraint associated with a 
solution to the homogeneous part of (I6.64p has been dealt with by introducing a global 
gauge fixing condition, so that no first-class constraint is set to zero strongly when we 
impose the second-class constraints to zero locally, and hence the Dirac bracket will be 
well defined. The Dirac bracket can be defined in the standard way. Then we can set the 
local second-class constraints IT- 7 , ^ and Tl^j to zero. Let us list the variables that 
are regarded as dependent variables: 

• The momenta p|- canonically conjugate to £* J can be set to zero. However, they did 
not appear in the Hamiltonian in the first place. 



The traceless component of £ y can be solved from \1/ 



i.r 



• \I/ can be regarded to constrain the metric gij. For instance we can regard that ^ 
fixes the conformal factor of the metric g^. Note that we cannot use \l/ to solve £ 
because \l/ is independent of £ lJ . 

• The constraints IP and tf/j can be regarded to define the traceless variables Qj and 
p^ in terms of independent variables. 

Now the remaining independent canonical variables are 11 variables in and p l \ the trace 
component £, and the trace components ( and p^. Surprisingly we could not remove £ 
from the set of independent variables even though the variables have a purely auxiliary 
role in the action (14.21) . 

Recall that in GR the canonical variables g^ and p^ are restricted by the Hamiltonian 
constraint and the momentum constraint so that we can freely specify Cauchy data for the 
trace and transverse-traceless parts of the momentum p %: > and the also spatial metric g^ 
up to a conformal factor 4^-42]. Then the longitudinal (vector) part of the momentum 



p 1 ^ and the conformal factor of metric are fixed by the momentum constraint and the 
Hamiltonian constraint, respectively. Fixing the trace p of the momentum fixes the slicing 
of space-time into spacelike hypersurfaces 4(J, |43j, |44j ; the scalar quantity g~ x ^ 2 p measures 



the rate of contraction of local three- volume elements with respect to local proper time. 
This leaves us the two well-known gravitational degrees of freedom which can be described 
by a pair of symmetric transverse traceless tensors. 

The remaining degeneracy left in the system is due to the time-dependent spatial 
diffeomorphism generated by the momentum constraint and the invariance under time 
reparameterization. Imposing the momentum constraint removes 3 degrees of freedom. 
Fixing the time reparameterization symmetry does not affect local dynamics. Thus the 
spatial metric tensor together with £ carries 3 physical degrees of freedom. The fact that 
the local Hamiltonian constraint is absent in the given theory implies that there should 
exist an additional scalar mode as in projectable version of HL gravity. But recall that 
here the projectability of N appeared because we chose to work with a foliation adapted 
to the field 4> so that the spatial hypersurfaces are those of constant 0. By solving the 



22 



equation of motion (12. 3p we have effectively chosen the direction of the time coordinate. 
Using an arbitrary spacelike foliation would not imply the projectability of N. Then of 
course there would appear higher derivatives of <fi in the action. In that case the extra 
mode would presumably be associated with itself. In addition there exists another extra 
degree of freedom in the trace component ( of the tensor field Qj. This extra mode is 
clearly associated with the second order time derivatives present in the action (13.29P . 
In the gauge N — 1, N l — 0, the Hamiltonian in reduced phase space is 

H = %= [ d 3 xU , (6.70) 



where 7i is expressed in terms of the independent variables in g^, p l i , £ and £, p^. 
However, we shall soon see that the Hamiltonian can be written without dependence on 
£, due to the constraints. 

The fact that a given higher derivative theory of gravity possesses extra degrees of 
freedom is always alarming, because such modes are often ghosts or otherwise pathologi- 
cal. Therefore it should be checked whether the extra modes present in this new version 
of CRG are ghosts, particularly the propagating mode ( of the higher-order variable. We 
can see that the kinetic part of the Hamiltonian contains terms involving the momentum 
which can attain arbitrarily negative values on the constraint surface. Note that this 
is not quite as evident as it would seem at first sight because of the existence of compli- 
cated constraints. Indeed a better view of this point is achieved by taking out a linear 
combination of the constraints from the Hamiltonian: 



i 



^ 2 p ij F^ kl V kl + _ 2n 2 F i r: 1 kl V kl g m g^F~ 1 nop V^ 



IK \ _L / ' 1 



C J l2n 2 a 2 U£ 



K 2 ( 

_| p-lpkl tm Jnp-1 pop 

8U 2 i) i) ± mnopS ^ 

+ 9aU ^ U; 3 C C (6.71) 

4 / \ 2 2 / 



( R 



2 - 2aU,C 3 D k D k Q lJ 



where we denote 

V kl =P M + 1 (( k m ( lm - ^CnnC") J ■ (6-72) 

This suggests that there exists a degree of freedom that carries negative energy. As is 
usual in higher time derivative theories, energy of the higher-order degree of freedom has 
opposite sign compared to the corresponding lower order degree of freedom. Here these 
gravitational degrees of freedom evidently interact with each other. We emphasize that 
the Hamiltonian is not the energy of the gravitational system, but rather a constraint that 
vanishes everywhere on the constraint surface. The ADM energy of an asymptotically flat 
gravitational system can be obtained by studying the surface terms of the Hamiltonian 
with appropriate boundary conditions. The problem is not that the Hamiltonian or the 
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total energy would attain arbitrarily negative values. The real source of the problem 
is that any state can decay into compensating positive and negative energy excitations. 
Even "empty space" will decay into a tempest of positive and negative energy excitations. 
This makes the theory unstable. The only way this could be avoided are the constraints. 
Unfortunately the Hamiltonian constraint is global and therefore it does not affect local 
physics. We suspect the momentum constraints does not protect the stability either. 

We note that the Hamiltonian formulation of CRG possesses a singularity at Qj = 0. 
In the subspace Qj = of phase space, all the momenta p 1 ^ are primary constraints and 
therefore there is no degree of freedom. Thus there is a kind of discontinuity between the 
solutions discussed in Sec. [5] and the ones with nonvanishing Qj. 

Some characteristics of the CRG action (I4.2p can be demonstrated by the following 
simple Lagrangian with two degrees of freedom q a (t) (a = 1, 2) 

L{q\ q a ) = \A ab q a q b - V(q a ) , A = "fj , (6.73) 

where a is a nonvanishing constant. The canonically conjugated momenta are defined as 

dL 

Pa = W = A ab q b . (6.74) 
When q 2 ^ 0, A has an inverse and we obtain the Hamiltonian 

^•>=^-K^) 2+W) - (675) 

In the limit q 2 — > 0, the kinetic part of the Hamiltonian diverges. Indeed when q 2 = 0, 
the conjugated momentum p2 = and the model reduces to a single degree of freedom 
q 1 . This is analogous to the case Qj = of CRG. When q 2 ^ 0, equations of motion are 

■ i P2 . dV 

aq z oqi 

q 2 = — - P2 p 2 = PlP2 - P2 - — (6 76) 

aq 2 (aq 2 ) 2 ' a(q 2 ) 2 a 2 (q 2 ) 3 dq2 

The Hamiltonian is not bounded below even when we consider solutions where q 2 retains 
its sign under time evolution. The kinetic matrix A has both a positive and a negative 
eigenvalue 

A± = \ (l ± v 1 + 4 (^ 2 ) 2 ) • ( 6 - 77 ) 
Thus we can write the Lagrangian as 

L = ±(Q+) 2 + ^(Q-) 2 -V(Q±) (6.78) 

± 



where the variables Q are defined by 



i 



A± „.±-a ,.+ / A± _ \ 2 I A- 



Q -<?' "=((5#+ 2 J (679) 

Thus the system has a positive energy and a negative energy degree of freedom. If this 
model were a continuum field theory it would certainly be unstable, assuming the two 
modes interact with each other. 
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We expect the CRG actions corresponding to higher values of the critical exponent 
z > 3 to exhibit similar problems as the case z = 3. As was already noted in Sec. [3j the 
number of time derivatives present in the ADM representation of the action grows with z. 
For example for z = 4, 5, 6 we are most likely to recover the same problems but in an even 
more complicated form than in the case z = 3. CRG actions corresponding to sufficiently 
high z will necessarily be unstable, once the order of time derivatives is greater than the 
number of available constraints. 

In more general perspective, we conjecture that generally covariant higher derivative 
theories of gravity which involve spontaneous (constraint induced) Lorentz and/or diffeo- 
morphism symmetry breaking will in general share the ghost problem with CRG. The 
only way this could be avoided is that either the higher time derivatives totally disappear 
(cancel out) from the ADM representation of the given action or the constraints available 
after symmetry breaking conspire to protect the stability of the higher-order degrees of 
freedom. The former way out is unlikely to be realized in d = 4 and higher spacetime 
dimensions because at least 2(d — 1) derivatives are required in the invariants present 
in the generally covariant action of the theory in order to achieve multiplicative renor- 
malizability. The latter way out is also hard to achieve especially when diffeomorphism 
invariance is broken and hence the constraint content is no longer similar to GR. Still it 
is not entirely impossible that such a theory could be conceived when our understanding 
on this kind of theories grows. 

7 Alternative variables for Hamiltonian formulation 

Instead of ( 14. ip alternative variables could be used to absorb the second order time deriva- 
tives in the action (I3.29P in order to construct Hamiltonian formulation of new CRG. Here 
we briefly consider two more choices of higher-order variables. 

7.1 Scalar variable 

Since the second order time derivatives in the action ( I3.29P can be written 



the following scalar variable can be used to absorb the second order time derivatives 



2 (/.»•' + KK*) d t fas + KKij) = <h (/»"' + KK*i) [I!,, + KKij)] 

+ AN (R* + KIT*) K lk g kl (R u + KK tj ) 
+ 4 (tfi + KK ij ) D {i N k) g kl (ify + KK tj ) 



(7.1) 



C = (R ij + KK ij ) [II 



(7.2) 



Now the action can be replaced by 



/ 



d 4 Xy/gN 



KjKi — K 2 + R 
2^2 




( 



4 (R ij + KK ij ) K lk K M {Ru + KK h ) 



+ a£ [C - {R ij + KK*) (Rij + KKij)] , (7.3) 
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where the scalar field £ is a Lagrange multiplier. This action has the disadvantage that it 
involves the extrinsic curvature to the sixth power. As a result the momenta canonically 
conjugate to the spatial metric will involve extrinsic curvature to the fifth power. Obtain- 
ing the Hamiltonian via Legendre transformation thus becomes a considerable problem. 
On the other hand the number of second-class constraints would be much lower than in 
the formalism of Sec. |6j 

However, the given form of the action is suitable for illustrating the question whether 
the theory contains ghosts or not. Let us consider the following expression that contains 
the only time derivative of the second order variable ( in the action 



*-5<*<k:;)( 



(7.4) 



Introducing two scalar fields 0i and 02 as 



j>i = ^=(K-Q, 2 = ^(^ + C), (7.5) 
enables us to diagonalize of this kinetic term. We obtain 

aObtfC = - #) ■ (7-6) 

Depending on the sign of the coupling constant a it seems that either <pi or 2 has a 
negative kinetic term which is an indication that given theory contains a ghost mode. 
But this is just the diagonalization of the cross-term between K and (. We also have to 
substitute 

K=^=(J> 1 + 2 ) , Ka = K l3 + M= (fa + 2 ) (7.7) 

into the hole kinetic part of the action (17.31) . One obtains that the action contains <\>\ and 
02 up to sixth power, including many cross-terms between <pi and 02- The sign of the 
kinetic terms with 0i and 02 to powers higher than second order is — sgn(a). In particular 
the terms that contains the highest powers are contained in the term 

-^(0 1 + 2 ) 6 . (T.8) 

Thus when a > the kinetic part does not possesses a stable minimum, while for a < 
it does. However, answering the question of ghosts may be more involved due to the 
presence of velocities to powers higher than two. Therefore the quadratic action ( 14. 2 p 
remains a safer alternative. 



7.2 Extrinsic curvature as an additional variable 

Yet another possible choice is to use the extrinsic curvature Kij as an additional variable 

&j = kj . (7.9) 

Then all the time derivatives in the action (13.291) will be absorbed into Qj and Qj. The 
time derivative of the trace of the extrinsic curvature can written 

K = g'Kij - 2iVC% - 2CW (t N 3) . (7.10) 
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The time derivative of the spatial Ricci tensor Rij can also be written in terms of Qj by 
using the relation 

g ij = 2NC ij + 2D {i N j) . (7.11) 

Then the action ( I3.29P can be replaced by one that depends on Qj, Qj and gij. Still the 
resulting action is more complicated than the one based on the choice (14. ip which is used 
in the Hamiltonian formalism developed in Sec. |6j 



8 Conclusion 

In this pap er we have analyzed a new version of the so called covariant renormalizable 
gravity 17J as an example of higher derivative theory of gravity which involves sponta- 
neous Lorentz and/or diffeomorphism symmetry breaking. It was earlier shown that this 
theory possesses the correct number of degrees of freedom when fluctuations around the 
flat background are analyzed and also that these fluctuations have modified dispersion 
relations so that the theory is power-counting renormalizable [TtJ. This is a remarkable 
feature, because the theory contains higher time and space derivatives which could imply 
potential problems with ghosts. For that reason it would certainly be nice to analyze 
the fully nonlinear theory with manifest diffeomorphism invariance. However, due to the 
complexity of the action we restricted our analysis to the case of the theory where the 
equation of motion (12. 3p is solved for the scalar field 0. Using a foliation of spacetime 
defined by <fi we obtained the ADM representation of the action of CRG and its Hamilto- 
nian structure. It turned out that the lapse function N has to depend on time only. As a 
result the theory obeys the projectability condition on N and lacks a local Hamiltonian 
constraint, which suggests that it is not possible to eliminate the additional scalar mode. 
Furthermore the fact that the CRG action involves higher time derivatives implies that 
it contains another additional degree of freedom, which can be seen when we introduce 
auxiliary fields. In summary, the theory contains 4 gravitational degrees of freedom. 3 
of the modes are carried by the spatial metric in conjunction with an auxiliary field £. 
Another extra scalar mode is carried by the trace component ( of the higher-order variable 
(14. ip which contains the second-order time derivatives present in the ADM representa- 
tion of the action (13.291) . We argued that the theory contains a degree of freedom that 
carries negative energy, which will destabilize the theory due to interactions with positive 
energy degrees of freedom. Apparently the available constraints are unable to prevent 
this pathology. As a result even "empty space" can decay into compensating positive and 
negative energy excitations. Then we conjectured that generally covariant higher deriva- 
tive theories of gravity which involve spontaneous (constraint induced) Lorentz and/or 
diffeomorphism symmetry breaking will in general share this problem with CRG. This 
is motivated by the difficulties in cancelling time derivatives from covariant higher-order 
theory and the reduced structure of constraints due to symmetry breaking. 

We should note that the secondary constraints of the theory are very complicated. 
Especially the second-class secondary constraints turned out inconveniently complicated. 
As a result the Hamiltonian has a very complicated form. This will make the canonical 
quantization of the theory a very hard task. Covariant quantization would likely be a 
feasible approach. 

Very recently we have begin to suspect that in this type of theory it might be better to 
treat the normal n M as a genuine dynamical variable with constraints restricting it to unit 
norm and zero vorticity. This way a more a general treatment could perhaps be achieved. 
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In covariant form the action would be defined as 



S 3 = / fxyf^ 



(4)D 

<g — - a(R al3 + KK afi - D a a p ) 

x (nW^ u + V^V M ) {R aP + KK aP - D a a p ) 

+ A (n^n" + 1) + B^T^ + M^B^B^j , (8.1) 

where the the vorticity for n M is 

J> = SV.V^] (8.2) 

and A, and M MJ , pcr are Lagrange multiplier fields. Variations of the action with respect 
to the Lagrange multipliers yields 

n/ = -l, 5^ = 0, .7> = 0. (8.3) 

The normal would associated with the scalar field by 

n M = -AV M (8.4) 

and choosing it to be the time = t (ADM gauge choice) yields the ADM formulation of 



the theory. This approach would be similar to the covariant formulation of HL gravity (45 
(also see [22, 13]), but with a much more complicated action: higher-order time derivatives 
and several extra kinetic terms etc. Indeed, the only conceivable advantage of CRG over 
HL gravity is the spontaneous breaking of Lorentz invariance in the high energy regime. 
The idea of achieving renormalizability of gravity via spontaneous symmetry breaking 
is certainly appealing, but we doubt whether it is worth the price of accepting such 
a complicated Hamiltonian structure, let alone a pathological extra degree of freedom 
which is simply unacceptable. 

Acknowledgements The support of the Academy of Finland under the Project Nos. 
136539 and 140886 is gratefully acknowledged. The work of J.K. was supported by the 
Grant agency of the Czech republic under the grant P201/12/G028. M.O. is supported 
by the Jenny and Antti Wihuri Foundation. 

A Appendix 
A.l Symmetrization 

Symmetrization and antisymmetrization of tensor indices is denoted by parentheses and 
brackets, respectively, e.g. 

B(a/3) = 7^ + Bp a ) , B[ a/3 ] = - (B a/3 — Bp a ) . 

We may also use the following notation if it is more convenient 

B aP + (a<+P) = 2B (a$) , B aj3 - (a fi) = 2B [a/3] , 
as can be the case in long expressions with many terms. 
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A. 2 Decomposition of the Riemann tensor 

Decomposition of the Riemann tensor of spacetime into components tangent and normal 
to the spatial hypersurfaces S t is given by the following identities. 

1. Gauss relation 

g\g\g p S^R\ pa = R\ aP + K\K 6fi - K\K aS (A.l) 

2. Codazzi relation 

s>v^y%v = 2d w k \] ( a -2) 

3. Ricci equation 

^X2> ff(4) ^V = K ^ K % + -ft D « D P N - ^NnK a p (A.3) 
These imply for example the decomposition of the scalar curvature of spacetime 

(4) i? = R + KijK ij -K 2 + 2V P (n»K) - —D i D { N , (A.4) 

where 

2 



2V M (n»K) = 2K 2 + -(k- N%K^j 



A.3 Decomposition of covariant derivatives of tensors tangent 
to spatial hypersurfaces 

Here we shall obtain space-time decompositions for first and second-order covariant deriva- 
tives of tensors tangent to spatial hypersurfaces E t . 

In the following calculations we will frequently need the covariant derivative of the 
orthogonal projector onto the spatial hypersurface 

V p g\ = (#/ - n p a u ) n a + n v (K^ - n p a a ) . (A.5) 

Consider a symmetric rank 2 covariant tensor field A a p that is tangent to E t , such as 
the two tensor fields inside the parentheses in (I3.24p . Such a tensor field is invariant under 
the orthogonal projector, A a $ = g^ a g v pA^, since n a A a/ 3 = n a A/3 a = 0. Let us decompose 
the covariant derivative of A a p: 

VpAap = 5y v (g p a g a pA pa ) = (g\ - n%) g" a g%V v A p<J + V M (g p a g%) V 
= g\g p a g a p VvA pa - n^g p a g^V n A p(7 + (V^g" a g% + g p a V^g^) A pa 

= DpAap - n M \JyC Nn A a p - 2K {a v Ap ) A + 2 (K* - n p a v ) n [a A p)v (A.6) 
= D^Aap - n M (—C Nn A a p - 2K {a u A l3)u + 2n {a Ap )u a v J + 2K^n {a A p)u 



29 



The component of V M A Q/3 that is fully tangent to S t is simply g l ' p g p a g 

For completeness let us list the rest 7 components of the decomposition of V^A^: 





V (7 T~~7 A 


1 r a 


-n Q 


g\n p g%V v A pa 


= -K/Ayp , 




9\9 P a n°V v A pa 


— — K v A 




n v n p g^V u A pa 


= —a u A u/3 , 




n u g p a n a V u A pa 


n u A 

Lb .n-OLV 1 




g\n p rfV v A pa 


= o, 


-n p n a rtp 


n v n p n a V u A p(T 


= 0, 



where in the left-hand side column we present the part of each component that are normal 
to S 4 in the decomposition (IA.6I) . 

For future needs we generalize Eqs. (13.181) . (13.191) and (IA.6[) for a generic covariant 
tensor field B ai ... Qk that is tangent to S t : 

^C Nn B av .. ah = V n B av .. ak +(K a / - n ai a p ) B n + - ■ - n ak a") B ttl .„„ , (A.7) 

l -X Nn B ai ... ak = g^ ai ■ ■■g p \ k V n B, 1 ..., k + K a /B n + ■■■ + K a /B ai ...„ , (A.8) 



N 



1 



V ^B ai ... ak D ^B ai ... ak n p | jy£>Nn,Ba.i—a.k B u ... ak ■ ■ ■ B ai ... u 



+ n a u B + h n a v B +n K V B -\ \- n K V B 

(A.9) 

It is worth to note that due to the implication a M = of the condition A" = N(t) [see 
(I3.22p ]. in ( 1A.7j) V n B ai ... ak is already tangent to E t , and hence the projection of (1A.7j) in 



f]A.8|) is trivial. However, as we did in (" 13 . 1 8[) and (IA.6j) . we continue to write a p explicitly 
in the decompositions of the covariant derivatives, in order to retain generality in this 
point. 

Then consider the second order covariant derivative of A a p. We decompose the second 
order covariant derivative of A a p similarly as the first-order covariant derivative in (1A.6|) : 



V,V u A aP = 5 p 5:V p V a {g x a g T p A XT ) 

= [g% - n p n p ) {g\ - n°n u ) g x a g^V p V a A Xr 

+ 2V (H (g\g T p ) V ]v) A XT + V P V U {g\g T p ) A Xr 

= g%g%g\g T ^ P ^aA XT - n p n p g\g x a g T ^ p V a A XT 

- n u g p ^g x a g T p V p V CT A XT + n p n u n p n a g\g T p V p V a A Xr 

+ 2V (H (g p a g%) V H A pa + V P V U {g p a g a p) A pu . 



(A.10) 



However, now we still have to decompose all the terms in flA.lOj) . In order to obtain the 



component of V '^V 1 v A a p that is fully tangent to S 4 we use (I3.25j) to calculate 
D p D v A afi = g p ^ v g x a g^V p D a A XT = g p ,g\g x a g T ^ p (g 4 '.g* \g\V 

= g p »g%g\g T ^ p (jV/ t ) V^A m + g p ,g a v g x a g T ^ p V a A M 
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which together with (1A.6|) finally gives 

g p p g a u g x a g T ^ p V a A Xr = D p D u A aP - K„ v (^C Nn A af3 - 2K {a p A p) )j + 2K^ al K/A pW . 

(A.11) 

Next consider the — n M component of the decomposition of V ^ v A a $ [the second term in 
(1A.10I) ]. It can be decomposed by using f !3.25|) to write 

V n D v Aap = V n (g%g\g T p V a A XT ) = V n (g%g x a g T fi ) V a A Xr + n p g%g x a g T p V p V a A Xr , 

(A.12) 

where the last term is what we are looking to decompose. Decomposition of the left-hand 
side of flA.12j) is given by flA.7j) 

V n D u A a)3 = —L Nn D v A a p - (K/ - n u a p ) D p A a/3 



(K a p - n a a p ) D v A pP - - n^a") D u A ap . 



Decomposition of the first term in the right-hand side of ( 1 A. 12p is given by (IA.6[) and 
taking the derivative of the projectors, 

Vn {g%g X a g T p) V ct A At = n v a p D p A aj3 + a u (^C Nn A aP - 2K {a p Ap )p j 

+ 2n {a \D u A Wp a p - 2K v p a^ a A f})p . 
Thus we obtain the — n p component of the decomposition (lA.lOp of V p V u A a/3 

n p g\g\g T p V p V a A XT = ±-C Nn D u A aP - K u p D p A a/3 - 2K (a p D v A Wp 

/ i \ (A.13) 

- a u I —C Nn A af3 - 2K {a p A l3)p J + 2K l/ p a {a A /3)p . 

The — n u component of the decomposition of V U V v A a p [the third term in flA.101) ] can be 
written 

g p p rfg\g^V p V a A XT = n p g\g\g T p V P V a A Xr + g p ^g\g^[V P) V a }A Xr . (A.14) 

In the second term we use the Ricci identity and the fact that A a p is tangent to S t to 
obtain 

g P Xg\g^[V p , V a ]Axr = ~g p p n°g\g^ ( ( %VV + (4) < p(7 ^) 

= -g p p rfg x a g r fi g^ {^% Xpa A^ T + ( % Tp(T ^) 

= -g» P n°g% [g x a {i) R p ^ x A\ + g T ^R p aVT Aj) 

= -2D [pl K, {a] A p p-2D [pl K pm A a p 
= 2D [a \K^A P m - 2D p K p{a A p p) 

where we have used the Codazzi relation flA.21) in the second last equality. Thus the —n v 
component of the decomposition flA.101) of "SI v A a p is 

g p ,n°g\g^V p V a A XT = ^C Nn D p A aP - K/D p A aP - 2K {a p D p A mp 

- a p (Jj£ Nn A a p - 2K {a p A p) )j + 2K/a (a Ap )p ( A - 15 ) 
+ 2D (a \K^A P m - 2D p K„ (a A P p ) . 
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Next consider the +n u n„ component of the decomposition of V P V u A a f3 [the fourth term 
in (jA.lOp ]. Let us write 

Vn (g\g T p V n A XT ) = V n (rfg\g T p ) V a A XT + nV^ V p V CT A Xr , (A. 16) 

where the last term is what we are looking to decompose. The left-hand side can alter- 
natively be written by using (IA.7j) and flA.8j) 



Vn (g x a g T pV n A Xr ) = V n \J^£ Nn A aj3 - 2K {a p A f})p ^j 



N 



-X 



Nn 



—C Nn A a p - 2K {a p Ap )p 



- (K a » - n a a p ) (^L Nn A pP - 2^^^ 
~ ( K p P ~ n P aP ) (jf£NnA ap - K { jA p) A . 

We also need to compute 

V n {n CT g x a g T p ) V a A Xr = a?D p A afi + n a a p \j^C Nn A pP - 2K {p °A P) A 

+ ripa p (j^C Nn A ap - K {a a A p) ^j - 2a (a Ap )p a p 
Thus we obtain the +n p n u component of the decomposition flA.101) of V„Vt,A a/ 3 
n p n a g x a g T /3 V p V a A Xr = j^C-Nn [j^C Nn A a A - ^£ Nn K {a P Ap )p 

~ 4K {Ql P ^£NnA W)p + 2K (a p K^A pa + 2K (a p K/A W)a (A ' 17) 

- a p D p A afi + 2a {a Ap )p a p . 

The rest 12 components of the decomposition of V u Vj,A a/3 are given in the last two terms 
of Eq. (lA.lOj) . Let us calculate them 



2V (H (g p a g%) V H A 



pa flat 



K/D v A pP - n p ( a p D„A pP + K u p -C Nn A 



pS 



- K/K p °A a0 - K u p K B a A pa ) + n u n„ ( a p -C Nn A 



1 



N 



1 pP 



- a p K°A a p - a p K°A pa ) + v 



+ [{KJ - n p a p ) (K va - n u a a ) + (//<-> v)] A p0 + (a (3) , 

(A.18) 

V„V„ {g p a g%) A pa = V^ u g p a A p0 + V p g p a V u g^A p(7 + (a^/3) 



n. 



n,. 



DpK up - K pu a p - n p ( —£ N nK„ n - 2K va K° n - a v a 



N 



Nn^up 



n p n v 



N 



iu {D p a p - K/K ap ) + 

+ [(K^ - n p a p ) (K ua - n u a a ) + (// o u)] A pfi 
+ n a n p (K / - n^a p ) (K* - n v a a ) A pa + (a 0) 



A p r _ 



(A.19) 
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Taking the sum of (1A.18j) and (1A.19|) gives 

2V M {g p a g%)V W) A pa + V M V, (g p a g° p ) A 



a p D v A oB + K„ p —£ Nn A oB + 1-XnuK up A p p 



1 pP 



v N 



N 



3K/K/A a g - K/K/A p(T - a u a p A 



ip/3 



- n u (^a p D p A pP + D p a p A p p + K/^C Nn A pB 



(A.20) 



2K^K°A aB - K/K/A pa 

1 



+ n p n u 2 [ a p —£ Nn A pj3 -r —^ Nn a p iM - z,w i\ p s± aj3 



--C-Nna p A p p - 2a p K p °A a 
- a'KfA^ + n a n p (K/ - n p a p ) {K v a - n y a°) A pa 

Now we can read the rest of the components of the decomposition of V p V v A a g from 



-n a g p p g a u n x g T pV p V a A XT = -2K { fD u) A pB - D p K up A p p + K pu a p A 



pP 



-np g p p g\g x a n T \7 p V a A XT = -2K (fM p D v) A ap - D p K up A a p + K pu A ap a p 



+n p n a n p g r7 u n x g T /3 V p V a A XT 
+n p n p n p g%g x a n T V p V a Axr 
+n u n a g p p n a n x g T B V p V a Ax T 



= -K/^£ Nn A pB - ^£ Nn K vp A p p - a p D„A 
+ 3K v p K°A aB + K/K^A p(7 + a v a p A pB , 

= ~~ K u p —£.NnA ap — —CN n K up A a p — a p D u A, 
+ 3K v p K p ° A aa + K v p K a a A pa + a u a p A ap , 

= -K l p —£ Nn A pB - a p D p A pB - D p a p A p p 
+ 2K p K o °A a0 + K p Kf A 



v 1 y a.p 



+n v np g p ^n a g x a n T V p V a A X r = -K fi p —£ Nn A ap - a p D p A ap - D p a p A 

+ 2K, p K n a A aa + K P K° A 



a 1 '■per j 



+n a n B g p n\n x n T V p V (T A XT = 2K P K/A 



flV V'- * P * AT 



per j 



n p n a n x g T B V p V a A X r = -2a p -j-£ Nn A pB - ^C Nn a p A p ^ 



+ Aa p K n a A aB + 2a p K°A 



-n p n v n B 

-n p n a n B 
-n v n a n g 
+n p n u n a n g 



g ■"■pa i 

n p n a g x a n T V p V a A XT = -2a p ^£ Nn A ap - ^£ Nn a p A a p 



+ Aa p K p a A a(T + 2a a K a p A pa , 



n p g%n*n T V p V a Axr = 2a p K/A pa , 
j'/n^VA = 2K lx p a a A pa , 
n p n a n x n T V p V a A X r = 2a p a a A pa . 
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Finally we have obtained every component of the decomposition of V^V^A^. 

One thing remains to be uncovered in the space-time decomposition. The Lie derivative 
(1A.7|) contains a time derivative of B ai ... ak that needs to be uncovered. For this purpose 
we assume a coordinate system on S 4 and use the connection coefficients explicitly. We 
denote the connection coefficients of V and D with 



and 



-pi 

1 jk 



-g im (djg m k + dkgjm — d m gjk) 



respectively. Using ( 13 .3p and ( 13 .4p we obtain 
(4) r° = 1 (N + N l d t N + KijrfNl 

(4) rg i = ^(^ + ^ i ) , 



(4) F _ j_ K 



(4) r* = g ij (Nd 3 N + N 3 - DjN k N f 



N' 

— [N + N J djN + K jk N J X 1 



(4) r 



n , - (Mk,, + DjN k ) 



~~N 



(djN + N k K kj ) 



(4)pi 



N*K. 



jk 



jk 



jk 



N 



(A.21) 
(A.22) 

(A.23) 
(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 



Then we obtain the covariant derivative in (1A.7|) by using the connection coefficients. For 
our purposes it is sufficient to consider only the spatial components of (I A .71) because the 
action (14 .2p involves only them (I4.3p . We obtain 



V n 5 n ... 



TV 



(Q B ■ - Wt v .b • Wr. R. 1 



1 

N 



B 



N J D 3 B n ... lk 



(NKX + D h N j ) B 
(NK ik j + D ik W)B iv . 

Thus we can write the spatial components of (1A.7I) as 



(A.29) 



N 



I'h n ~ N j DjB iv .. ih - D h N^Bj... ik 



N 



Di^B^ 



(A.30) 



where Cn denotes the Lie derivative along the shift vector N l in H t 
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